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We present a numerical model study of the zero-temperature infrared optical properties of 
(III,Mn)V diluted magnetic semiconductors. Our calculations demonstrate the importance of treat- 
ing disorder and interaction effects simultaneously in modelling these materials. We find that the 
conductivity has no clear Drude peak, that it has a broadened inter-band peak near 220 meV, and 
that oscillator weight is shifted to higher frequencies by stronger disorder. These results are in good 
qualitative agreement with recent thin film absorption measurements. We use our numerical find- 
ings to discuss the use of f-sum rules evaluated by integrating optical absorption data for accurate 
carrier-density estimates. 

PACS numbers: PACS numbers: 73.20.Dx, 73.20.Mf 



I. INTRODUCTION 

Recent studies of (III,Mn)V diluted magnetic semicon- 
ductors (DMSs), motivated by the discovery of carrier- 
mediated ferromagnetism jl| in these materials, have un- 
covered a wealth of interesting phenomenology involving 
an interplay between collective electronic effects, broken 
symmetries, interactions, and disorder. The richness of 
the transport and optical properties of (III,Mn)V DMSs 
jl], ^| arises from their strong valence-band spin-orbit cou- 
pling and from the sensitivity of their magnetic state to 
growth and annealing conditions, doping, and external 
fields H, Q]. In situ tunability of material properties by 
applying gate voltages ||, by illumination ||], or by ap- 
plying external magnetic fields, suggests many opportu- 
nities for exploring new physics or building new devices 
using these materials. There is considerable evidence 
that the properties of these systems can be understood 
using a simplified effective model [Q in which the low 
energy degrees of freedom are S = 5/2 Mn 2+ local mo- 
ments, and holes in the host-semiconductor valence band. 
The Hamiltonian of this effective model should in general 
include terms that represent scalar and exchange inter- 
actions between the Mn ions and the valence band holes, 
and Coulomb interactions among valence band holes and 
between the holes and the Mn ions. The simplest the- 
oretical approach to these models, one that appears to 
be relatively successful in explaining many bulk proper- 
ties of strongly metallic DMS ferromagnets, either ignores 
disorder or treats it perturbatively and approximates ex- 
change interactions using a mean-field theory |t], ||, [| [l(J 
p5[ . We have recently presented a theory of the intra- 
band optical conductivity of (III,Mn)V ferromagnets that 
uses such an approach In the present paper we re- 
port a finite-size numerical study of the optical conduc- 
tivity of (Ga,Mn)As in which disorder is treated exactly, 
enabling us to address weakly metallic or insulating sys- 
tems and to understand some of the qualitative strong- 
scattering effects that are omitted in the simpler theory. 



Our interest in this property of DMS ferromagnets is 
motivated by a series of recent optical absorption exper- 
iments |l^, [l^, [l4|, [l5|]. The experiments, performed in 
thin film geometries in both metallic and insulating sam- 
ples (as judged by the temperature dependence of the re- 
sistivity below T c ), exhibit several common features: (a) 
non-Drude behavior in which the conductivity increases 
with increasing frequency in the interval between meV 
and 220 meV, (b) a broad absorption peak near 220 — 260 
meV that becomes stronger as the samples are cooled, 
and (c) a broad featureless absorption between the peak 
energy and the effective band gap energy which tends to 
increase at the higher frequencies. The 220 meV peak 
can be attributed to inter- valence-band transitions [ pd| , 
to transitions between the semiconductor valence band 
states and a Mn induced impurity band, or to a com- 
bination of these contributions. Dynamical mean-ficld- 
theory studies jl6), for a single-band model that neglect 
the spin-orbit coupling and the heavy-light degeneracy 
of a III-V semiconductor valence band, demonstrate that 
non-Drude impurity-band related peaks in the frequency- 
dependent conductivity occur in DMS ferromagnet mod- 
els when the strength of exchange interaction coupling is 
comparable to the valence band width. The conductiv- 
ities predicted by this model are, however, inconsistent 
with experiment in temperature trend and appear not 
to be in the same regime as experimental samples. On 
the other hand we Jll[ demonstrated in earlier work that 
peaks in the conductivity can also occur in the weak- 
coupling regime, provided that the realistic multi-band 
character of the valence bands is acknowledged. In this 
approach the 220 meV peak is due to heavy-hole to light- 
hole inter-valence-band transitions and the gradual in- 
crease in absorption observed in the experiments |l3| at 
frequencies above 500 meV is ascribed to transitions to 
the bands split-off by strong spin-orbit interactions. Al- 
though this theory could account for many overall fea- 
tures of the measured optical conductivity in metallic 
samples, one prominent feature was at odds with ex- 
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perimental data, namely the relative magnitudes of the 
lu —>■ conductivity and the 220 meV conductivity peak. 
The numerical calculations reported on in this paper sug- 
gest that this discrepancy is due to multiple-scattering 
effects that suppress the low-frequency conductivity. 

We restrict our attention in this paper to the T = 
limit, which allows us to neglect scattering of valence 
band quasiparticles off thermal fluctuations in the Mn 
ion spin orientations. We assume a simple ground state 
in which the Mn spins have parallel spin orientations p6| . 
Sources of disorder that are known to have some impor- 
tance in (ni,Mn)V ferromagnets include randomness in 
the placement of Mn ions on cation sites of the host lat- 
tice, random placement of Mn ions at interstitial sites, 
and randomly located antisite defects in which group V 
elements are placed at group III element lattice sites. In 
the low-energy effective model, Mn ions at cation sites 
have charge Q — — e while interstitial Mn and antisite 
defects have Q — 2e. Coulomb and exchange interac- 
tions with the Mn ions have comparable importance |Tl[] 
for the scattering of valence band holes and must both 
be included in the theory. In the approach of Ref. [If] 
these sources of disorder scattering were treated pertur- 
batively using a Born approximation, evaluating the va- 
lence band quasiparticle lifetimes using Fermi's Golden 
rule |T^| . In this paper we find the T = quasiparticle 
states for finite-size systems by solving the Schrodinger 
equation. Unlike the Born approximation theory, this ap- 
proach does not assume metallic behavior and remains 
valid at low Mn concentrations when jl8| Mn acceptor 
impurity bands emerge. We find that corrections to the 
Born approximation theory are important at a quantita- 
tive level even for strongly metallic high Mn concentra- 
tion samples, and that they explain the property that the 
low-frequency and dc conductivities are suppressed com- 
pared to their values in the Born approximation theory. 
We also find that the f-sum rule, obtained by integrating 
the conductivity over frequency, differs from its Born ap- 
proximation value by less than 10% for typical metallic 
carrier densities. 

This paper is organized as follows. Section II defines 
the model Hamiltonian we use to describe the valence 
band system in a (III,Mn)V ferromagnet which adds ex- 
change and Coulomb interactions to the six band k • p 
envelope function Hamiltonian of the host (III,V) semi- 
conductor. Section III describes technical issues associ- 
ated with the evaluation of Kubo linear-response-theory 
formula for the finite-size system optical conductivity. In 
Sec. IV we discuss the numerical results we have ob- 
tained for different Mn and carrier density regimes, which 
are in qualitative agreement with experiment and discuss 
the degree to which aspects of the extreme impurity-band 
jl6| and Golden-rule gl], [l7| limits are reflected in the ex- 
act results. We also emphasize the feasibility of using the 
f-sum rule to extract reasonably accurate measurements 
of the carrier concentration Jlf| , circumventing the inher- 
ent uncertainties in Hall measurements in systems with 
large anomalous Hall effect coefficients. 



II. MODEL HAMILTONIAN 

We approximate the host semiconductor valence band 
by its six-band k • p Kohn-Luttinger model. This simpli- 
fication takes advantage of the small number of holes per 
atom in these ferromagnets. The band electrons interact 
via Coulomb and exchange interactions with randomly 
distributed Mn spins and with randomly distributed an- 
tisite defects via Coulomb interactions alone; we do not 
explicitly account for the possible role of Mn interstitials 
|I9| . The requirement of overall charge neutrality im- 
plies that the density of valence band holes (p) is related 
to the density of Mn (nMn) and the density of antisites 
(tias) by p = n Mn - 2nAs- In MBE grown (III,Mn)V 
ferromagnets, it is often a challenge to accurately deter- 
mine the density of holes, and even the density of Mn 
ions that have been incorporated into the ferromagnetic 
state, although progress is being made. It appears at 
present that any annealing procedure |IJ |2^] that leads 
to higher carrier densities, inevitably leads to higher fer- 
romagnetic transition temperatures, and to much higher 
conductivities For Gao.95Mno.05As for example, it 
is now possible to prepare samples with T = dc con- 
ductivities that are close to an order of magnitude larger 
than in the samples for which ac conductivity measure- 
ments have so far been carried out. Although we believe 
that much of the variabilty in physical properties can be 
accounted for in terms of variations in the fraction of 
Mn ions that are incorporated and, more importantly, 
variations in the carrier density, it appears likely that a 
part of this variation is still associated with sources of 
unintended disorder that are not present in our model. 
We expect that as unintended disorder is reduced, the 
properties of these ferromagnets will be more completely 
determined by the two parameters of our model calcula- 
tions, p and tiM n [0 ■ 

Below we will evaluate the frequency-dependent con- 
ductivity by using the Kubo formula description for the 
linear response of quasiparticles to an electromagnetic 
field. Coulomb interactions play an essential role in de- 
termining the quasiparticle states and interactions among 
the holes cannot be ignored. In this paper we use a 
Hartree mean-field approximation for the quasiparticles, 
which is consistent with the standard bubble diagram 
approximation we use for the conductivity. 

The single particle part of the quasiparticle Hamilto- 
nian may be written as the sum of host crystal band and 
disorder pieces H° = H h + H D : The disorder Hamilto- 
nian has three contributions H D = R K -^- + tfMn-h + 
H As ~ h . H Mn ~ h describes the attractive Coulomb inter- 
action between the ionized Mn 2+ acceptor and a valence 
band hole, H K ex - describes the local kinetic-exchange in- 
teractions between valence band carriers and the Mn mo- 
ments, while H As ~ h describes the repulsive Coulomb in- 
teraction between a hole carrier and ionized anti-site de- 
fects. The single-particle envelope function model Hamil- 
tonian then reads: 
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where J(r) = ( J pd ) / {{2Tia 2 f/ 2 )e- r2 ' 2a o , and J* = 
(s" x , s y , s z ) where s x ^ VyZ are the 6x6 matrices which de- 
scribe hole spins in the representation of the six band 
envelope-function Hamiltonian, and I is a 6 x 6 unit ma- 
trix. In this Hamiltonian, / labels Mn sites, K labels As 
sites, and Si stands for a Mn spin with quantum number 
S = 5/2. The Mn and As positions are denoted by Ri 
and Rk. The term proportional to Vq is a central cell 
correction [El) to the Coulomb attraction between holes 
and Mn acceptors that is known to be required to recover 
experimental values for the isolated acceptor energy. The 
quantities Nmu and Na s are the number of Mn and As 
atoms, which we distributed randomly inside cubic simu- 
lation cells of side L. Charge neutrality in our simulation 
cell means that the number of quasiparticle levels occu- 
pied by holes is Nh = Nmu — 2Na s - The values chosen 
for the phenomenological parameters that appear in this 
equation are J v d — 0.05eV nm 3 , ao — 0.3 nm, e = 10.9, 
V = 2.5eV and r = 0.259 nm. 

The host band part of the Hamiltonian is described 
via the six band Kohn-Luttinger model, with wavevec- 
tors k in the envelope function replaced by k a — iV a 
for a = x, y, z to account for the inhomogeneity induced 
by disorder. Choosing the angular momentum quan- 
tization direction to be along z-axis, and ordering the 
j = 3/2 and j = 1/2 basis functions according to the list 
(-3/2, 1/2, -1/2, -3/2; 1/2, -1/2), the Luttinger Hamil- 
tonian H L has the form IlOI: 
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In the matrix (0) we have highlighted the j = 3/2 sector. 
The Kohn-Luttinger eigenenergies are measured down 
from the top of the valence band, i.e. they are hole en- 
ergies. For completeness we list the expressions which 
define the quantities that appear in H L : 
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with 7i = 6.98, 72 = 2.06 and 73 = 2.93. 

The Hartree potential due to interactions among holes 
must also be included in the single-particle Hamiltonian 
for quasiparticle states since it captures the screening of 
Coulomb interactions with impurities due to the build up 
or depletion of charge. The total one-particle Hamilto- 
nian is H = H° + V H , where 



kj,k'j' 



EE 



Aire 2 

i? 2 '' l/r 



,k'+p+kf ■ 



(4) 



and the density matrix 
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with f a = 1(0) for an occupied (unoccupied) state. Here 
^£■(0 = e%k Tu k ( 7? )' wnere u l~-(f) are the Bloch func- 



tions. The Hartree potential must be determined by solv- 
ing self-consistcntly for the density matrix. We have di- 
agonalized this Hamiltonian using a plane- wave represen- 
tation for the envelope functions of H° and applying peri- 
odic boundary conditions in a cube of size L, which limits 
wavevectors to the discrete set k — (n x , n y , n z )(2n/L), 
where n x ,n yi n z are integers. For different system sizes 
we choose different maximum values of n XjVlZ , denoted by 
n m , so that the maximum wavevector k max is held fixed. 
For the calculations presented here k max — tt nm -1 . 



III. KUBO FORMULA FOR THE OPTICAL 
CONDUCTIVITY 

We evaluate the real part of the frequency dependent 
conductivity from the Kubo formula expression which re- 
lates it to the quasiparticle eigenvectors \a > and eigen- 
values E a . 
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Matrix elements of the momentum operator, p a = jV a 
with a = x,y, z, are given by 

< a\p a \? > = EE4f c % < *&ip«i**v > 



3>3 k,k' 



/-^ kj 



kj kj' h Qk a 



(7) 



At zero temperature and positive u) the difference (f a — 
fp) restricts a to occupied quasiparticle states and (3 to 
empty quasiparticle states. For finite-size systems the 
real part of the conductivity is more reliably evaluated 
by broadening the 5 function (2^] 
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with the level broadening, 7 , being of the order of the 
mean level spacing ~ Ep/N^ ~ 1/L 3 . In the thin film 
limit the real part of the conductivity is proportional to 
the absorption coefficient in transmission experiments, 
a.{u>) — with n = ^/cq being the zero frequency 

limit index of refraction of the thin film material and eo 
the dimensionless dielectric constant. This formula for 
a(uj) is most accurate in the far infrared regime or in the 
thin film limit since at higher frequencies multiple scat- 
tering from the interfaces must be taken into account. 
This can be easily done using the standard electromag- 
netic description of light waves propagating through a 
medium, but the simple proportionality is lost. 



IV. RESULTS 



A. Optical conductivity of parabolic band with 
disorder 



It is useful for pedagogical purposes to first present 
frequency-dependent conductivities calculated with this 
approach for the case where the 6-band Luttinger model 
is replaced by a parabolic band model with a mass equal 
to the heavy hole mass of GaAs. In this way we can par- 
tially separate the effects of strong disorder on the optical 
conductivity from the complications associated with the 
heavy, light, and split-off valence bands. In the calcu- 
lation all Mn spins point along the z-axis so that spin 
up and down single-particle Hamiltonians decouple. For 
each spin the Hamiltonian is identical to that of Eq. 1 
except that H L is replaced by the usual kinetic term 
-fi 2 V 2 /2m/ l ft and the 6x6 identity matrix is replaced by 
unity. The following parameters were used in the calcu- 
lation of the conductivity: p=0.33 nm -3 , nMn = 1 nm -3 
(x w 4.5%), L = 8 nm, and n m = 4. The results pre- 
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FIG. 1: Absorption and conductivity of minority (dotted 
line) and majority (dashed line) carriers in a parabolic band 
approximation. Solid line is the total absorption. Here 
p = 0.33 nm -3 , n Mn = 1 nm" 3 (a; w 4.5%). 



sented below were obtained by averaging over Nd = 10 
realizations of the disorder potential. 

Fig. showss conductivities and absorption coefficients 
we have evaluated for minority and majority carriers. 
Even for this simple band model the conductivity dis- 
plays clear non-Drude behavior, increasing as a function 
of frequency at low frequencies. We believe that the non- 
Drude behavior below 200 meV reflects multiple scatter- 
ing effects that enhance high-angle scattering and would 
lead to localization if the disorder were stronger. These 
effects are explicitly neglected in the Boltzmann trans- 
port theory which leads to the Drude formula. The f- 
sum rule for the conductivity calculated by integrating 
these numerical results over frequency deviates from the 
parabolic band model value by less than 10%, demon- 
strating that the errors induced by constructing a smooth 
conductivity curve using broadened S- functions are un- 
der reasonable control. The f-sum rule is more accurately 
satisfied for larger system sizes, for which we are able 
to choose smaller values of the broadening parameter 7. 
The non-Drude behavior we find here does capture one 
aspect of the experimental results. However the over- 
all conductivity magnitude of this model is much smaller 
than in experiment. Evidently the mixing of heavy and 
light hole bands leads to quasiparticles with larger veloc- 
ities. 



B. Optical conductivity of Luttinger-Kohn 
Hamiltonian with disorder 



Mixing of heavy- and light-hole bands and its interplay 
with magnetic order and scattering disorder can be in- 
vestigated by solving the full Hamiltonian H — H° + V H . 
This approach treats the band structure realistically and 
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FIG. 2: Absorption and conductivity of a metallic sample 
computed using a Luttinger Kohn Hamiltonian with disorder. 
Nominal parameters same as in Fig. pi. 



FIG. 3: Absorption and conductivity of a disordered sample 
computed using a Luttinger Kohn Hamiltonian. The parame- 
ters are p = 0.2 nm -3 , nun = 0.88 nm -3 , L = 6 nm, No = 4 
and n m — 3. 



the randomly located charged impurities without ap- 
proximation. Fig. |^ shows the frequency-dependent 
conductivity we have evaluated for p — 0.33 nm -3 , 
n-Mn = 1 nm -3 , L — 6 nm, — 4 and n m = 3. 

For hui < 200 meV we observe that the conductivity in- 
creases with frequency, unlike the Drude formula case. 
The increase is actually sharper than in the parabolic 
band model, a property that we believe reflects the ad- 
ditive contributions of non-Drude intra-band absorption 
and broadened inter-valence-band absorption. There are 
significant differences between these results and those of 
the Born approximation (see Fig. 1 of Sinova et al. 1 1 ) 
calculated for the same model parameters. In our cal- 
culation the maximum value of the absorption is about 
90000 cm -1 while in the calculation of Sinova et al. it 
is about 110000 cm -1 . We believe that this reduction is 
due to the enhanced backscattering effects captured by 
the exact calculation and also due to a transfer of spectral 
weight out of the intra-band contribution. The peak near 
hw ~ 0.9 eV in these calculations is due to heavy-hole to 
split-off band transitions and is more pronounced here 
than in the Born approximation calculations. The pres- 
ence of this peak at rather large frequencies emphasizes 
one weakness of the present calculation, namely that it 
neglects valence to conduction band transitions. It will 
be interesting to see whether or not future experiments 
confirm our prediction of a peak in the absorption due to 
these transitions that lies at high frequencies, but still 
clearly below the onset of valence-to-conduction band 
transitions. 

The results shown in Fig. || are for parameters that 
are typical for high T c (Ga,Mn)As ferromagnets. The 
dc limits of the conductivities that we calculate are sub- 
stantially higher than those of the samples for which cr(w) 
values have been published at present, but comparable to 
or smaller than those measured for the highest conductiv- 



ity (Ga,Mn)As ferromagnets currently available. Fig. || 
displays the absorption curves calculated for a more dis- 
ordered system |2q| . The parameters are p = 0.2 nm -3 , 
n-Mn = 0.88 nm~^ L = 6 nm, Nd = 4 and n m = 3. The 
conductivity is substantially suppressed at all frequen- 
cies, but most strongly at low frequencies. It again has 
non-Drude behavior at low- frequencies, increasing with 
increasing frequency over the interval between meV and 
220 meV. The broad heavy-to-light absorption peak near 
220 meV is still clear but the heavy-hole to split-off band 
peak is less visible, perhaps explaining the fact that it 
has not yet been identified experimentally. Hirakawa et 
al. report that the maximum value for absorption in the 
more disordered sample they study is about 20000 cm -1 , 
which has the same order of magnitude as the value cal- 
culated for the parameters of Figure [| 

These results are qualitatively consistent with overall 
features of the measured optical conductivity, including 
the overall mag nitude £! H [l§. The accuracy of our 
numerical results is limited somewhat by the wavevector 
cutoffs that we must invoke in order to make the numer- 
ical calculations manageable. This is especially trouble- 
some in estimating the importance of heavy-hole to split- 
off band transitions. The gradual increase in absorption 
beyond 500 meV seen in experiment may reflect these 
transitions. We expect that allowing more wavevectors 
in our calculations would shift oscillator strength toward 
higher energies and produce a broader peak, in better 
agreement with experiments. (Increasing n m from our 
current value, 3, by just one changes the dimension of the 
Hamiltonian matrix from 2058 to 4374, and the resulting 
self-consistent computation becomes quite substantial.) 
On the other hand, our investigation of parabolic band 
case, where n m can be increased to 4, suggests that sys- 
tem size dependence of low energy non-Drude part of 
absorption is negligible for n m bigger than 3. 
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V. CONCLUSIONS 

We have presented a non-perturbative self-consistent 
study of the optical properties of diluted (III,Mn) V DMS 
which treats both disorder and interactions on equal foot- 
ing. Our calculations for the diagonal ac-conductivity are 
in good qualitative or semi-quantitative agreement with 
optical absorption experiments in thin film geometries. 
Several of the non-Drude behaviors observed in the ex- 
periments are naturally accounted for in this approach 
which is able to capture multiple scattering effects that 
are beyond the scope of perturbative Born approximation 
approaches plf . Our calculations show that spin-orbit 
coupling of valence band holes must be included in the 
calculation in order to account for the qualitative features 
of the experiments. It is possible that features associated 
with heavy-hole to split-off band absorption will emerge 
more clearly near 0.9 eV when ct(oj) measurements are 
performed on the most metallic (Ga,Mn)As ferromagnets 
currently available. 

We have used these calculations to test the efficacy 
of the f-sum rule we proposed in earlier work for carrier 
density measurements. In applying the f-sum rule to the 
present system a complication arises, compared to micro- 
scopic and parabolic band cases, in that the integrated 
conductivity depends on the character of the occupied 
quasiparticle states. In previous work pj we estimated 



that for carrier densities in the range that is interest- 
ing for (Ga,Mn)As ferromagnets, the effective mass that 
should be used in the f-sum rule is ~ 0.24m e , intermedi- 
ate between heavy and light-hole masses. This estimate 
was based on calculations of the occupied quasiparticle 
states that treated disorder perturbatively in the Born 
approximation. We have found by integrating the opti- 
cal conductivities evaluated here over frequency that this 
optical effective mass estimate is reliable at the ~ 10% 
level, when disorder is treated exactly. Uncertainties in- 
troduced by our finite momentum cut-off do not allow us 
to estimate the value or the variability of this mass more 
accurately. Nevertheless, this level of accuracy should be 
comparable with what is possible from Hall effect mea- 
surements which are complicated by the strong anoma- 
lous Hall effect in these ferromagnets. 
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